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Estimation  of  Time  to  Maximum  Supercooling 
During  Dynamic  Frazil  Ice  Formation 


STEVEN  F.  DALY  AND  KATHLEEN  D.  AXELSON 


INTRODUCTION 

Frazil  ice  is  the  major  form  of  ice  in  northern  rivers  and  streams.  It  has  a  significant 
negative  impact  on  winter  navigation,  hydropower  production  and  the  use  of  water  in¬ 
takes,  and  contributes  to  ice  jam  flooding.  Three  rather  general  stages  of  frazil  evolution  can 
be  identified.  The  first  is  the  dynamic  nonequilibrium  stage,  characterized  by  supercooled 
water,  turbulent  flow,  rapid  growth  of  disk-shaped  crystals  and  the  creation  of  new  crystals 
by  secondary  nucleation.  The  second  stage  follows  the  first  in  time  and  is  the  evolution  and 
transport  stage.  The  third  stage  is  characterized  by  stationary  floating  ice  covers. 

The  first  stage  has  been  comprehensively  described  and  reproduced  in  the  laboratory, 
and  can  be  quantitatively  discussed  through  the  use  of  a  number  continuity  equation  and 
a  heat  balance  equation  (Daly  1984).  These  equations  are  dimensionally  incompatible  and 
strongly  nonlinear.  Simplified  forms  of  these  equations  are  used  with  success  to  predict  the 
steady-state  performance  of  industrial  crystallizers.  However,  the  dynamic  stage  of  frazil 
ice  is  inherently  a  nonsteady  process  and  has  resisted  analysis. 

Merrier  (1984)  developed  a  global  model  to  describe  reactive  transport  of  dissolved  and 
suspended  solids  in  lakes,  rivers  and  coastal  areas.  In  an  extension  of  the  work  of  Daly 
(1984),  Merrier  applied  his  model  to  frazil  ice.  In  addition,  Mercier  suggested  certain  sim¬ 
plifications  that  could  be  applied  in  limited  situations.  To  explain  the  sensitivity  of  the 
global  model  to  variations  in  the  system  parameters,  for  example,  Mercier  determined  an 
analytical  expression  for  the  time  required  for  a  turbulent  mass  of  water  undergoing  con¬ 
stant  heat  loss  and  a  constant  seeding  rate  to  reach  the  point  of  maximum  supercooling.  The 
expression  retains  the  basic  physics  and  produces  an  analytical  result  amenable  to  simple 
numerical  solution.  In  this  paper  we  will  describe  the  development  of  this  analytical  ex¬ 
pression,  report  the  numerical  solution,  review  laboratory  experiments  in  which  the  time 
to  reach  maximum  supercooling  was  measured,  explore  the  influence  of  the  system  para¬ 
meters,  apply  the  expression  to  the  laboratory  results,  and  discuss  the  implications  for  frazil 
growth  in  rivers  and  streams. 


THEORY 

The  analytical  expression  proposed  by  Mercier  assumes  a  homogeneous,  isotropic 
water  body.  Initial  frazil  production  is  assumed  to  be  caused  by  seed  crystals,  whose  size 
is  near  the  critical  radius,  inuuduceu  into  Ihe  water  at  a  constant  seeding  rate.  The  valu_  of 
the  critical  radius  can  be  calculated  using  thermodynamic  principles  (Forest  and  Sharma 
1987)  and  is  estimated  to  be  about  0.4  gm  for  the  range  of  supercooling  found  in  rivers  and 
streams  (Mercier  1984).  The  expression  assumes  a  completely  mixed  system,  with  uniform 


temperature  and  seed  crystals  quickly  transported  into  suspension  throughout  the  water 
column.  Additional  crystals  are  produced  by  secondary  nucleation.  It  is  assumed  that 
secondary  nucleation  is  dominated  by  new  crystals  created  through  collisions  of  large 
crystals  by  the  action  of  turbulent  shear.  The  size  of  the  new  crystals  is  assumed  to  be  close 
to  the  critical  radius. 

The  thermal  energy  balance  is  determined  in  the  following  manner.  The  water  is 
supercooled  by  a  constant  heat  loss  rate  per  unit  volume.  Above  the  equilibrium  tempera¬ 
ture,  any  seed  crystals  introduced  will  melt.  Below  this  temperature,  seed  crystals  can  grow 
and  produce  new  crystals  by  secondary  nucleation.  It  is  assumed  that  as  the  ice  crystals 
grow,  the  latent  heat  of  fusion  released  is  negligible  before  the  point  of  maximum 
supercooling  is  reached.  At  the  point  of  maximum  supercooling,  the  rate  of  latent  heat  of 
fusion  released  equals  the  heat  loss  rate.  After  this  point,  the  water  temperature  rises  as  the 
frazil  ice  growth  continues.  The  point  of  the  maximum  increase  in  temperature  corresponds 
to  the  maximum  rate  of  frazil  growth.  Figure  1  illustrates  this  process. 

A  convenient  parameter  to  characterize  frazil  ice  growth  is  the  time  until  maximum 
supercooling  fs.  For  particles  less  than  100  am  in  size  (which  is  a  good  assumption  for  frazil 
crystals  before  the  maximum  supercooling),  frazil  ice  crystal  growth  prior  to  maximum 
supercooling  may  be  described  by  the  following  equation 


dr  _  hT 
dt  pi  L 


where  r 
t 
h 
T 

Pi 

L 


crystal  radius 
time 

heat  transfer  coefficient 
supercooled  water  temperature 
density  of  ice 
latent  heat  of  fusion. 


(1) 


A  good  approximation  for  the  heat  transfer  coefficient  (Daly  1984)  is 

h  -  L  (2) 

where  k  is  the  thermal  conductivity  of  water.  Assuming  that  the  latent  heat  released  by  the 
growing  frazil  is  negligible  up  until  the  time  of  maximum  supercooling,  Mercier  described 
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Figure  2 .  Typical  curve  of  water  temperature  versus  time  during  frazil  ice 
growth. 
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the  temperature  of  the  supercooled  water  by  the  following  equation 


(3) 


where  Q  =  heat  loss  rate  per  unit  volume 
Cp  =  heat  capacity  of  the  water 
p  =  density  of  water. 

Substituting  eq  3  into  eq  1,  and  assuming  that  the  critical  radius  r.is  approximately  0,  we 
find  the  crystal  radius  r  at  any  time  t  using  the  equation 


lit)  = 


\pp/  LCp ) 


(4) 


where  X  is  the  coefficient  of  crystal  growth.  1  he  general  equation  for  secondary  nucleation 
is  (Daly  1984) 


It  =  ZEt  (5) 

where/T  =  number  of  secondary  nuclei  produced  per  unit  time  per  unit  volume 
Z  =  number  of  crystals  produced  per  unit  of  collision  energy 
Et  =  total  rate  of  energy  transfer  by  the  various  mechanisms  of  collision. 

The  mechanisms  can  include  collisions  caused  by  fluid  shear,  turbulent  shear  and  buoy¬ 
ancy,  and  collisions  with  boundaries.  For  the  present  case,  it  is  assumed  that  only  collisions 
between  two  crystals  are  important,  and  eq  5  may  be  written  in  the  form 


ItU)  =  Z 


dn,  r\)dn  dr\ 


(6) 


where  r.and  rare  the  radii  of  the  two  particles  colliding  and  c(r,r)  is  the  collision  function 
(Mercier  1984). 

Taking  secondary  nucleation  into  account,  Mercier  gives  a  simplified  particle  number 
density  distribution  as  follows 


rc<  r<  Xt 
r>  Xt 


(7) 


where  g(r,t)  is  the  crystal  number  density  distribution  and  Jv  is  the  crystal  seeding  rate  per 
volume  of  water.  To  approximate  the  number  of  crystals  produced  per  unit  collision 
energy,  it  is  assumed  that  the  primary  cause  of  collisions  is  turbulent  shear,  and  that 
collisions  of  the  largest  particles  (where  r  ~  \t)  will  be  the  most  important.  At  r  ~  Xt,  g(r,t) 
will  be  approximately  1JX.  Assuming  rc  ~  0,  this  approximation  yields  for  secondary 
nucleation 


g(r,0  =  Ir(t)  =  0.018Zp,(E-)3/2(^)2[l.67-0.451nAf](Af)10  (8) 

v  \  X  I 

where  e  is  the  turbulent  energy  dissipation  rate  and  v  is  the  kinematic  viscosity  of  water  (see 
Mercier  [1984]  for  details).  Noting  that  the  mass  of  ice  can  be  written 

M,(t. )  =  p,  nr3g(r,t)dr  (9) 
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and  given  the  approximation  of  the  number  density  distribution  of  eq  7,  after  substituting 
from  eq  8,  the  mass  of  ice  is 

Mj(f )  =  |  pi  X3  fv  f 4  +  (i  7t)2  (0.018)  Zp2g-)3/2/v2  Xu(4.3  xlff4)  t 14.  (10) 

At  the  time  of  maximum  supercooling, 

dMi  =  0_  (11) 

dt  L  ' 

Then,  by  finding  the  derivative  of  eq  10  and  solving  for  fs,  an  expression  for  the  time  until 
maximum  supercooling  is 

tid+Ql0)  =  (12a) 

Pi  IAJ/v 

where  £  is  a  coefficient  defined  as 

5  =  4.5x  KT4  piZ  (J-)3/2/vX8  (12b) 

and  X  is  defined  in  eq  4.  Equation  12  is  an  analytical  expression  for  the  time  to  maximum 
supercooling  in  terms  of  four  system  parameters:  heat  loss  rate,  crystal  seeding  rate  per  unit 
volume,  turbulent  energy  dissipation  rate  and  number  of  crystals  produced  per  unit 
collision  energy.  Two  of  the  parameters,  heat  loss  rate  and  turbulent  energy  dissipation 
rate,  can  be  measured  and  are  controlled  in  every  frazil  ice  experiment.  The  seeding  rate  is 
much  more  difficult  to  control  and  measure.  The  number  of  crystals  produced  per  unit 
energy  is  a  material  property  that  has  been  estimated  for  other  materials,  but  can  only  be 
estimated  by  numerical  experiment  for  frazil  ice. 

The  suitability  of  this  analytical  expression  may  be  verified  as  the  time  to  maximum 
supercooling  has  been  reported  for  several  different  experiments.  So,  we  solved  eq  12  and 
evaluated  the  relative  influence  of  the  parameters.  The  results  were  compared  to  the  experi¬ 
mental  data  as  shown  below. 


SENSITIVITY  OF  ANALYTICAL  SOLUTION  TO  SYSTEM  PARAMETERS 

Using  a  computer  program  that  we  developed,  we  solved  eq  12  for  ts.  Heat  loss  rate, 
crystal  seeding  rate  per  unit  volume,  number  of  crystals  produced  per  unit  collision  energy 
and  turbulent  energy  dissipation  rate  were  varied.  V alues  of  the  other  parameters  used  are 
shown  in  Table  1 .  To  verify  the  accuracy  of  the  program,  the  results  were  compared  to  those 
reported  by  Mercier.  The  comparison 
showed  that  the  program  performed 
satisfactorily. 

Following  verification,  we  evalu¬ 
ated  the  relative  influence  of  each  para¬ 
meter  on  the  time  to  maximum  super¬ 
cooling  by  varying  each  over  several 
orders  of  magnitude,  but  within  a  real¬ 
istic  range.  The  results  are  shown  in 
Figure  2  for  a  constant  heat  loss  rate, 
with  seeding  rate,  turbulent  energy 
dissipation  rate  and  number  of  crys- 


Table  1.  Values  of  parameters  used  to  determine 
sensitivity  of  the  analytical  solution  to  system 
parameters  (after  Batchelor  1967). 


Parameter _ Symbol _ Value 


Density  of  ice 

P, 

0.92  g/cm1 

1  g/cm1 

Density  of  water 

P 

Latent  heat  of  fusion 

L 

330 J/g 

Thermal  conductivity 

k 

5.6  xIO'J/cm  s  'C 

Specific  heat  capacity 

4.217  J/g  *C 

Kinematic  viscosity 

V 

1.787  x  10  :cmVs 

4 


c.  e  =  1  cm2 /s*. 


Figure  2.  Effects  of  varying  system  parameters  on  the  time  to  maximum 
supercooling  Cts>  for  a  heat  loss  rate  Q  of  10^1  /cm3 s  and  three  turbulent 
energy  dissipation  rates  e. 


tals  produced  per  unit  collision  energy  varying  as  indicated.  The  time  to  maximum 
supercooling  is  dominated  by  the  crystal  seeding  rate  Is  and  heat  loss  rate  Q.  The  number 
of  crystals  produced  per  unit  collision  energy  Z  was  less  of  a  factor  unless  the  crystal 
seeding  rate  was  low. 

The  number  of  crystals  produced  per  unit  collision  energy  has  a  greater  effect  on  time 
to  maximum  supercooling  at  higher  turbulent  energy  dissipation  rates  (high  e).  This  is 
reasonable  since  more  collisions  would  be  expected  at  higher  turbulence  and  hence 
secondary  nucleation  would  play  a  greater  role  in  frazil  ice  formation,  particularly  if  the 
initial  seeding  rate  is  small. 


EXPERIMENTAL  DATA 


Frazil  ice  experiments  have  been  conducted  under  a  variety  of  conditions,  and  time  to 
maximum  supercooling  has  been  reported  for  several  experiments.  In  each  experiment, 
although  the  heat  loss  rate  and  turbulent  energy  dissipation  are  set  by  the  experimenter,  the 
values  of  these  parameters  have  not  often  been  reported.  However,  based  on  the  informa¬ 
tion  provided,  it  is  possible  to  calculate  these  values  as  described  below. 

The  heat  loss  rate  may  be  determined  from  the  rate  of  temperature  decline  of  the 
turbulent  water  prior  to  the  formation  of  frazil.  To  determine  the  turbulent  energy 
dissipation  rate,  the  Reynolds  number  of  the  flow  in  the  experiment  is  computed  from  the 
average  velocity  and  a  characteristic  length.  For  pipe  flow,  the  diameter  of  the  pipe  is  the 
characteristic  length,  and  for  open  channel  flow,  the  hydraulic  radius  is  the  characteristic 
length.  The  friction  factor  is  determined  from  the  Reynolds  number  of  the  flow  using  an 
estimated  surface  roughness  (Chow  1959).  For  pipe  flow,  the  turbulent  energy  dissipation 
rate  may  be  found  using 


e 


_  r/2i 

2  d 


(13) 


where  /  =  friction  factor 

it  =  average  velocity 
d  =  pipe  diameter. 

For  open  channel  flow',  the  turbulent  energy  dissipation  rate  is  a  function  of  the  friction 
velocity  it.,  which  is  the  product  of  the  friction  factor  and  the  mean  velocity 


ii. 


(14) 


The  turbulent  energy  dissipation  rate  is  then 


£ 


(15) 


where  K  is  von  Karman's  constant,  generally  set  equal  to  0.4,  and  Rhis  the  hydraulic  radius. 

Table  2  contains  a  summary  of  data  from  the  frazil  ice  experiments  described  below’. 
Figure  3  depicts  the  range  of  values  for  heat  loss  rate  and  turbulent  energy  dissipation  rate 
found  from  the  experimental  data. 

Michel  (1963)  reportedly  carried  out  over  80  frazil  ice  growth  experiments  in  an  outdoor 
recirculating  flume  constructed  of  a  Plexiglas  channel,  30.5  cm  wide  by  30.5  cm  deep  by  6.7 
m  in  length,  connected  to  a  15-cm-diameter  Plexiglas  pipe.  A  20-cm-diameter  variable 
speed  pump  with  cast  iron  casing  provided  circulation.  Many  of  the  tests  were  done  at 
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Table  2.  Summary  of  experimental  data. 


Test  t  Q  f 


Exverimenter 

Description  of  experiment 

no. 

(min) 

(liem' s) 

(iin:i<‘) 

Michel  (1963) 

Outdoor  recirculating  flume. 

27 

3.3 

1.2x10' 

33 

Carstens  (1966) 

Recirculation  flume  in-10'C 

6A 

4 

1.19x10  1 

13 

coldroom;  propeller  in  tank 
in  -1()‘C  coldroom. 

5A 

8 

3.9xl()‘ 

Hanlev  and 

Cylindrical  tank  with  pad- 

_2’c 

38.1 

0.9  x  10' 

.3 

Michel  (1977) 

dies  in  coldroom.  Varied  cold- 

-5*C 

21.6 

1.9  x  10'1 

3 

room  temperature  and  vel- 

-1()‘C 

14.7 

3x  H)'1 

3 

ocitv. 

-20‘C 

8.1 

7.1  x  10-' 

3 

Tsang  and 

Warm  air  jacket  around 

A2  (48'<< ) 

118.8 

1.9  x  It)1 

2 

Hanlev  (1985) 

tank,  except  test  C,  which 

Bl-3  (23'7, ) 

43.6 

2  x  10  1 

was  in  recirculating  flume 

B2-3(U'.,7> 

36.4 

2  x  10' 

Varied  salinity.  Seeded 

B3— 4  (fresh) 

18.1 

1.9  x  10' 

-> 

with  shavings  of  ice.  Aver¬ 
age  coldroom  temp.  -1  C. 

C-4  (ocean) 

123.9 

1.1  x  10  ' 

0.7 

Mueller  (1978) 

Supercooled  water  first  in 

E08 

— 

— 

1375 

agitation  tank  in  warm  room. 

E09 

— 

— 

4Wv 

night,  and  air  temperature  ranged  from  -32  to  -6.7*C.  Water  was  pumped  into  the  flume 
and  then  circulated.  Water  temperature  was  measured  every  15  seconds  by  a  differential 
thermometer,  reportedly  accurate  to  ±  0.0025°C.  Data  were  reported  for  only  one  test,  in 
Michel's  Figure  2.  The  water  temperature  decline  rate  was  reported  and  we  calculated  the 
heat  loss  rate  from  this.  Assuming  that  the  experiment  was  dominated  by  the  turbulence 
that  occurred  immediately  downstream  of  the  pump,  we  calculated  the  turbulent  energy 
dissipation  rate  for  this  flow  area.  The  reported  average  flow  was  375  gal./ min  (0.025  nr1/ 
s),  yielding  a  velocity  in  the  pipe  of  73  cm/s,  a  Reynolds  number  of  8.3  x  104  and  a  friction 
factor  of  0.023.  The  turbulent  energy  dissipation  rate  was  then  found  using  eq  1 3.  Time  to 
maximum  supercooling  was  measured  from  the  figure. 

A  recirculating  oval  flume,  20  cm  wide  by  30  cm  deep  by  600  cm  long,  located  in  a  -10"C 
coldroom,  was  used  by  Carstens  (1966)  to  study  frazil  ice.  Water  depth  was  20  cm.  The 
bottom  and  sides  of  the  acrylic  flume  were  insulated,  and  cooling  of  the  water  was  aided 
by  a  fan  blowing  along  the  straight  portion  of  the  flume.  Circulation  in  the  flume  was 
induced  by  a  variable  speed  propeller.  Water  temperature  was  measured  by  a  hand-held 
mercury  thermometer  marked  to  0.01°C,  immersed  5  to  10  cm.  In  a  typical  test,  the  water 
temperature  decline  rate  was  measured  and  the  heat  loss  rate  was  then  determined .  For  the 
experiments  reported,  the  average  flow  velocity  was  50  cm  / s,  yielding  a  Reynolds  number 
of  1 .8  x  104  and  friction  factor  of  0.019.  We  calculated  the  turbulent  energy  dissipation  rate 
using  eq  15.  Unfortunately,  for  several  experiments  shown  in  Carsten's  Figure  5a,  the 
turbulent  energy  dissipation  ;ate  could  not  be  calculated  because  the  velocity  of  the  water 
was  not  reported.  However,  the  time  to  maximum  supercooling  could  be  determined  from 
this  figure  and  is  included  for  comparison. 

Hanley  and  Michel  (1977)  conducted  frazil  ice  experiments  in  a  stainless  steel  tank,  120 
cm  in  diameter  and  76  cm  deep,  at  air  temperatures  of -2,  -5,  -10  and  -20°C.  Water  depth 
in  the  tank  was  not  reported.  Paddles  attached  to  an  axis  and  located  25  cm  above  the 
bottom  of  the  tank  induced  rotational  velocity  in  the  tank.  For  these  experiments,  the 
coldroom  temperature  was  set,  and  then  the  paddles  started.  Water  temperature  was 
measured  by  a  thermistor  placed  just  below  the  water  surface.  We  calculated  the  heat  loss 
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Figure  3.  Estimated  values  of  turbulent  energy  dissipation  rate  and  heat  loss  rate 
for  experimental  data  and  three  rivers. 


rates  from  the  reported  water  temperature  decline  rates.  The  results  presented  by  Hanley 
and  Michel  were  reportedly  averaged  over  all  velocities;  the  minimum  velocity  at  which 
frazil  growth  was  noted  was  24  cm/s,  and  the  maximum  reported  velocity  was  73  cm/s. 
Assuming  a  water  depth  in  the  tank  of  70  cm  and  an  average  velocity  of  48.5  cm/s,  we 
calculated  the  Reynolds  number  at  8.8  x  104  and  estimated  the  friction  factor  to  be  0.0145. 
Equations  13  or  15  yielded  the  same  turbulent  energy  dissipation  rate.  Time  to  maximum 
supercooling  was  not  reported.  Total  time  of  supercooling  was  reported,  however,  as  was 
the  time  from  the  minimum  temperature  to  the  end  of  supercooling.  We  could  estimate 
time  to  maximum  supercooling  by  subtracting  the  latter  from  the  former. 

Frazil  ice  formation  in  fresh  water,  ocean  water  and  artificial  seawater  of  varying  salinity 
was  examined  by  Tsang  and  Hanley  (1985).  They  placed  a  rectangular  Plexiglas  tank,  38  cm 
long  by  25.5  cm  wide  by  15  cm  deep,  equipped  with  a  stirrer  at  one  end  to  provide 
turbulence,  in  a  -15’C  coldroom.  A  horizontal  plate,  0.75  times  the  length  of  the  tank,  was 
placed  at  mid-depth  to  produce  vertical  recirculation  currents,  and  a  jacket  of  air,  slightly 
above  0"C,  was  placed  around  the  bottom  and  sides  of  the  tank.  Cooling  was  provided  by 
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a  fan  that  pushed  air  across  the  surface  of  the  water.  Tsang  and  Hanley  measured  water 
temperature  with  a  thermometer  reportedly  calibrated  to  O.OOOl'C,  with  repeatability  to 
0.001  °C.  In  these  experiments,  the  water  was  first  cooled  to  the  previously  selected 
temperature  and  then  seeded  either  by  a  ball  of  ice  or  by  scrapings  from  ice.  We  picked  one 
sample  from  each  of  their  five  groups  of  experiments  for  analysis  here.  We  calculated  heat 
loss  rates  for  each  sample  as  we  did  for  Hanley  and  Michel.  All  but  the  group  C  experiments 
took  place  in  the  tank  described  above,  with  a  reported  average  velocity  of  15  cm/s.  Water 
depth  was  not  reported  and  we  assumed  it  to  he  1 5  cm.  For  these  experiments,  we  calculated 
the  Reynolds  number  to  be  3.9  x  103,  with  a  resulting  friction  factor  of  0.031 .  We  estimated 
turbulent  energy  dissipation  rates  using  eq  13.  Group  C  (ocean  water)  tests  took  place  in  a 
recirculating  flume,  15  cm  wide  by  13  cm  deep,  with  a  water  depth  of  1 1  cm,  and  average 
velocity  of  15  cm/s.  They  reported  the  Reynolds  number  as  8.54  x  10\  From  this,  we 
estimated  the  friction  factor  to  be  0.025.  The  turbulent  energy  dissipation  rate  was 
estimated  using  eq  15.  Time  to  maximum  supercooling  was  reported. 

Mueller  (1978)  studied  the  nucleation  process  of  frazil  ice  using  an  agitating  tank,  17.2 
cm  long  by  12  cm  wide  by  20  cm  deep,  in  a  coldroom  kept  slightly  above  0”C.  The  tank  was 
surrou  nded  by  a  jacket  through  which  a  coolant  was  circulated.  Agitation  was  provided  by 
a  grid  submerged  in  the  tank.  In  these  experiments,  the  water  was  supercooled  to  the  de¬ 
sired  degree,  the  coolant  circulation  was  stopped,  agitation  was  begun  and  the  supercooled 
water  seeded.  Water  depth  ranged  from  20  to  20.5  cm.  He  recorded  temperature  using 
linear  type  thermistors  that  were  accurate  to  a  reported  +0.02’C  in  air  and  ±0.002°C  in  water. 
Turbulent  energy  dissipation  rates  were  reported.  These  experiments  were  not  designed  to 
allow  the  time  to  maximum  supercooling  to  be  determined  as  the  water  was  not  seeded 
until  the  maximum  supercooling  was  achieved .  W e  report  the  turbulent  energy  d issipa  tion 
rates  only  to  provide  a  comparison  with  the  other  experiments. 


APPLICATION  TO  EXPERIMENTAL  DATA 

We  applied  the  analytical  expression  to  the  experimental  data  cited  above,  assuming  a 
constant  number  of  crystals  produced  per  unit  energy  of  4  x  10'°  nuclei/erg  as  suggested 
by  Mercier.  Given  the  experimental  values  of  heat  loss  rate  and  turbulent  energy  dissipa¬ 
tion  rate  for  each  experiment  (Table  2),  we  optimized  the  seeding  rate  to  produce  the 
reported  times  to  maximum  supercooling.  The  results  are  shown  in  Figure  4  and  listed  in 
Table  3.  The  optimized  values  of  lv,  which  range  from  10~'  to  10~*  crystals/cm3  s,  appear 
to  be  reasonable  based  on  the  experimental  conditions. 

A  more  consistent  means  of  comparing  the  optimized  seeding  rates  is  to  compare  the 
results  as  a  seeding  rate  per  unit  surface  area,  rather  than  a  seeding  rate  per  unit  volume. 
The  seeding  rate  per  unit  surface  area  can  be  simply  calculated,  given  the  design  of  the 
experimental  apparatus.  Michel's  experiment  had  the  highest  opti-  rized  surface  seeding 
rate  (I  )  of  2.7  cr  Hs/ cnr  s.  Although  the  rate  of  seeding  was  not  measured  or  considered, 
it  is  not  diffir.:’  .  imagine  that  conducting  an  experiment  outdoors,  as  these  were,  would 
be  a  situatic  hich  a  large  number  of  seed  crystals  could  be  present.  The  experiments 
of  Carstens  condo  "  ’.ed  in  a  large  coldroom,  had  a  comparably  high  value  for  the  optimized 
surface  seeding..  .  .11  crystal/  cm2  s.  The  presence  of  Carstens'  recirculating  oval  flume, 
driven  by  <  mopuler,  in  this  coldroom  could  have  contributed  to  the  production  of  seed 
crystals.  The  optimized  seeding  rates  for  the  experiments  of  Hanley  and  Michel  are  less 
than  those  of  either  Michel  or  Carstens.  However,  the  results  for  these  experiments  show 
a  consistent  increase  in  seeding  rate  with  a  decrease  in  the  coldroom  air  temperature.  While 
this  type  of  relationship  has  not  been  measured  or  reported  elsewhere,  such  an  inverse 
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Figure  4.  Calculated  time  versus  reported  time  for  experimental 
data  using  optimized  seeding  rate  I... 


Table  3.  Comparison  of  experimental  and  computed  results. 


Exfvri- 


Source 

Q 

<l/cm' s) 

£ 

mental 

p 

(min) 

Model 

t, 

(min) 

ls  b* 

/crystals  \  An/ttals 
V  cm's  /  V  cm’s 

Michel  (1963) 

1.2  x  10  1 

33 

3.3 

3.0 

10'1  2.7 

Carstens  (1966) 

1.19x  10  ’ 

13 

4 

3.7 

5  x  10':  1.0 

Hanlev  and  Michel  (1977) 

-2'C 

9x  10-1 

3 

38.1 

40.5 

10'5  7  x  l()-‘ 

-5'C 

1.9  x  10  J 

3 

21.6 

21.6 

10-1  7x10-' 

-10’C 

3x  10  1 

3 

14.7 

13.0 

10  '  7  x  1()‘: 

-20*C 

7.1  x  10-* 

3 

8.1 

7.5 

5x10''  3.5x10'' 

Tsang  and  Hanlev  (1985) 

48',', 

1.9  x  10  ’ 

2 

118 

<  10" 

23r<< 

2x  10  ’ 

2 

43.6 

<  10" 

1 1  %c 

2x  10  1 

2 

36.4 

<  10-6 

fresh 

1.9  x  10  1 

2 

18.1 

18.3 

5x10"  7.5x10' 

ocean 

1.1  x  10  ' 

0.7 

124 

<  10" 

*  Michel,  assumed  depth — 27  cm;  Carstens,  reported  depth — 20cm;  Hanley  and  Michel, 
assumed  depth — 70  cm;  Tsang  and  Hanlev  assumed  depth — 15  cm. 


relationship  between  seeding  rate  and  room  air  temperature  seems  likely.  Figure  5  is  a  plot 
of  the  coldroom  air  temperature  and  the  calculated  seeding  rate  for  each  experiment  for 
which  these  values  were  reported  or  could  be  calculated.  Air  temperature  for  the  experi¬ 
ment  given  in  Michel  was  not  reported;  therefore,  the  average  of  reported  air  temperatures 
for  all  tests  is  plotted. 

The  optimized  seeding  rates  for  the  experiments  conducted  by  Tsang  and  Hanley  are 
very  low,  less  than  10"* crystals/cm1  s.  We  do  not  know  how  much  the  design  of  their  ex¬ 
periment  contributed  to  these  optimized  seeding  rates.  The  experiments  were  conducted 
in  a  box  surrounded  on  all  sides,  except  for  the  open  top,  by  a  warm  air  jacket.  A  horizontal 
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Temperature  (°C) 


Figure  5.  Plot  of  optimized  seeding  rate  (I.)  versus  experimental  air 
temperatures. 

baffle  separated  a  top  layer  and  a  bottom  layer.  It  may  have  been  that  the  seed  crystals  were 
not  present  in  their  coldroom,  or  that  the  crystals  were  not  fully  mixed  in  their  tank,  or  that 
the  seed  crystals  melted  near  the  walls  of  the  tank. 

Unfortunately,  it  is  not  really  possible  to  know  how  realistic  these  optimized  values  are 
as  there  are  no  measured  values  reported  either  from  the  laboratory  or  the  field .  Osterkamp 
(1977)  reports  that  hexagonal  plate  ice  crystals  ranging  from  60  to  350  pm  in  size  were 
observed  on  the  air  above  a  supercooled  Alaskan  stream.  Their  concentrations  ranged  from 
6  to  6  x  104  crystals/ m\  It  is  possible  to  estimate  a  surface  seeding  rate  from  the 
concentration  of  crystals  suspended  in  the  air  CAas 

K  =  Ca",  16) 

where  »(is  the  terminal  fall  velocity  of  the  ice  crystals.  Estimating  a  terminal  fall  velocity 
of  the  largest  and  smallest  ice  crystal  as  17 cm/s  and  5  cm/s,  respectively  (Pruppacher  and 
Klett  1980),  yields  a  seeding  rate  per  unit  surface  area  of  3  x  10~?to  1.0  crystals/cnr  s.  This 
range  encompasses  the  optimized  surface  seeding  rates  (Table  3)  and  provides  a  certain 
degree  of  support  for  them. 


TIME  TO  MAXIMUM  SUPERCOOLING  IN  NATURAL  WATER  BODIES 

No  field  data  on  time  to  maximum  supercooling  exist  for  rivers.  In  fact,  the  temperatures 
of  supercooled  water  are  rarely  reported.  To  investigate  the  implications  of  the  analytical 
analysis  for  rivers  and  streams,  we  examined  three  frazil  producing  rivers:  the  Ottauque- 
chee  River  in  Vermont,  the  Susquahanna  River  in  Pennsylvania,  and  a  1 60-km  reach  of  the 
Salmon  River  near  Salmon,  Idaho.  Using  average  values  for  channel  slope,  width,  depth, 
discharge  and  Manning's  n,  we  calculated  the  average  velocity  for  each  river  assuming  uni¬ 
form  flow.  The  heat  loss  rate  was  calculated  assuming  an  average  heat  transfer  rate  of  20 
W/m2  ’C  for  the  Ottauquechee  and  Susquehanna  Rivers  and  10  W /nr  °C  for  the  Salmon 
River,  as  suggested  by  Zufelt.*  The  estimates  of  turbulent  energy  dissipation  rate  are 


*  Personal  communication  with  J.  Zufelt,  CRREL,  1988. 
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Table  4.  Analytical  solution  applied  to  three  frazil  producing  rivers, 
Z  =  4  X  1010  nuclei/erg,  Iy  in  crystals/cm3s. 


Time  (min) 


River 

Qtllcni' s) 

eh.Hr/s1) 

I..  =  0.00  7 

1,  =  o.ot 

1=0.1 

1,=  1 

Salmon 

1  x  10-' 

50 

13.8 

9.6 

6.3 

3.4 

Ottauquechee 

3.28  x  10-1 

44 

9.3 

6.5 

4.4 

2.7 

Susquehanna 

1.31  x  10-* 

3.5 

17.0 

11.6 

6.9 

3.3 

shown  with  the  experimental  data  in  Figure  3.  Using  the  computer  program,  we  calculated 
time  to  maximum  supercooling.  The  results,  shown  in  Table  4  for  various  seeding  rates, 
indicate  that  time  to  maximum  supercooling  in  rivers  is  quite  short,  on  theorder  of  minutes. 

Considering  the  magnitude  of  the  estimated  times,  field  measurement  of  time  to  max¬ 
imum  supercooling  will  be  quite  difficult.  For  example,  assuming  a  crystal  seeding  rate  per 
unit  volume  of  10-1  ciystals/ cm5  s,  the  time  to  maximum  supercooling  for  the  Salmon  River 
would  be  about  6  minutes.  At  an  average  winter  velocity  of  about  1 .73  m/s,  supercooling 
would  occur  within  a  reach  of  660  m.  This  length  is  about  0.4%  of  the  1 60-km  reach  of  river 
in  question.  This  may  have  implications  for  the  placement  of  ice  control  structures  that  rely 
on  adfreezing  of  frazil  ice. 

SUMMARY 

Mercier  (1984)  presented  a  complex  model  of  frazil  ice  growth.  Our  investigation  ex¬ 
amined  a  simplified  analytical  expression  that  described  the  time  until  maximum  super¬ 
cooling,  a  parameter  often  measured  in  laboratory  experiments,  in  terms  of  four  system 
parameters:  heat  loss  rate,  seeding  rate,  turbulent  energy  dissipation  rate  and  the  number 
of  crystals  produced  per  unit  of  collision  energy.  To  explore  the  influence  of  the  system 
parameters,  a  computer  program  to  solve  for  time  to  maximum  supercooling  was  devel¬ 
oped.  The  program  was  verified  using  Mercier's  data.  Using  a  fixed  value  of  the  number 
of  crystals  produced  per  unit  collision  energy  suggested  by  Mercier  (4  x  10ll)  nuclei/erg), 
and  the  experimental  values  of  heat  loss  rate  and  turbulent  energy  dissipation  rate,  we  em¬ 
ployed  the  program  to  optimize  the  seeding  rate  to  produce  the  reported  times  to  maximum 
supercooling.  This  process  resulted  in  times  close  to  those  reported  for  reasonable  values 
of  seeding  rate.  An  analysis  of  the  data  indicated  that  the  seeding  rate  and  heat  loss  rate  have 
the  most  influence  on  time  to  maximum  supercooling,  while  the  number  of  crystals  pro¬ 
duced  per  unit  energy  and  the  turbulent  energy  dissipation  rate  have  less  influence. 

We  found  an  inverse  relationship  between  the  coldroom  air  temperature  at  which  the 
experiment  was  conducted  and  the  seeding  rate.  Using  estimated  values  of  heat  loss  rate 
and  turbulent  energy  dissipation,  we  computed  a  range  of  values  for  time  to  maximum 
supercooling  for  the  Salmon,  Ottauqueechee  and  Susquehanna  Rivers.  The  results  indicate 
that  time  to  maximum  supercooling  in  rivers  and  streams  may  be  very  short,  and  that 
supercooled  regions  may  be  difficult  to  locate  and  measure. 
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